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STUDIES  OF  THE  EFFECT  OF 
DEPTH  OF  FO^JS  ON  SEISMIC  PULSES 

I .  Abstract 

This  report  summarizes  the  cor. elusions  reached  in  thirty 
other  reports  issued  in  connection  with  this  contract,  and 
lists  the  individual  reports  which  detail  how  each  of  the 
objectives  of  the  contract  was  reached. 

The  principal  method  used  in  an  attempt  to  measure  depth- 
of-focus  of  shallow  events  is  the  Watson-Merdler  inverse¬ 
filtering  technique.  It  is  shown  that  a  simplified  version 
of  this  method  which  neglects  noise  also  finds  ghosts  in  a 
seismogram,  as  does  cross-correlation.  The  three  methods 
frequently  but  not  always  find  the  same  ghosts. 

Two  important  factors  which  make  it  difficult  for  all 
three  methods  to  find  ghosts  are  that  short-period  seismometers 
do  not  record  the  main  part  of  the  energy  in  most  seismic 
pulses  and  the  ghosts  commonly  differ  greatly  from  their 
primaries  in  waveform  and  frequency  spectrum. 

II.  Objectives 

The  objectives  of  this  research  are  listed  below  to¬ 
gether  with  reference  to  the  scientific  and  technical 
reports  which  describe  each  part  of  the  work.  Inasmuch  as 
detailed  technical  reports  have  already  been  distributed  on 
most  aspects  of  this  research,  only  the  conclusions  will  be 


|JPP|rr"ll lllinnni — 1 - -vmmmpm  T~“  . . 


2 


summarized  here  except  for  the  work  of  the  last  three  months, 
which  has  not  been  thoroughly  reported  previously. 

Objective  1:  To  construct  seismic  scale-models  wherein 
the  effects  of  the  earth  on  seismic  pulses  can  be  simulated. 
The  model  facilities  are  described  in  Technical  Report  Nr. 

2  by  R.  L.  Rothman  (1964).  Figure  1  shows  typical  arrange¬ 
ments  of  some  of  the  equipment.  Scientific  Report  Nr.  2  by 
Y.  Nakamura  (1963A)  discusses  the  behavior  of  some  of  the 
mode]  materials.  Scientific  Report  Nr.  3  by  B.  F.  Howell, 

Jr.  (1966)  describes  a  digitizer  developped  for  use  in  this 
laboratory. 

Objective  2:  To  study  interference  patterns  in  model 
records  and  to  compare  them  with  patterns  observed  in  real 
events.  Technical  Report  Nr.  2  by  Rothman  (1964)  is  con¬ 
cerned  with  interference  in  the  model  records.  Technical 
Report  Nr.  4  by  G.  Leblanc  (1966)  discusses  the  spectra  of 
a  real  earthquake  in  detail.  This  problem  has  also  been 
discussed  by  Nakamura  (1963B)  in  Technical  Report  Nr.  1. 

Objective  3:  To  collect  seismograms  for  the  study  of 
depth  of  focus.  A  file  of  over  400  seismcgrams,  most  of 
them  in  digital  form  has  been  assembled  as  a  source  of 
research  material. 

Objective  4;  To  try  various  methods  of  analyzing  seis¬ 
mograms  in  an  attempt  to  measure  accurately  the  depth  of 
focus  of  the  causitive  event.  This  was  the  principal  ob¬ 
jective  of  the  project.  It  was  hoped  that  accurate  depth 


determination  would  assist  significantly  in  distinguishing 
natural  earthquakes  from  blasts.  A  method  of  inverse  filter¬ 
ing  called  the  Watson-Merdler  method  was  developped,  and  is 
described  in  Technical  Report  Nr.  3  oy  S.  C.  Merdler  (1964). 
The  original  method  was  modified  and  tested  on  blasts  as 
described  in  Technical  Report  Nr.  5  by  J.  L.  Lin  (1966)  and 
on  earthquakes  as  described  in  Technical  Report  Nr.  6  by 
Y.  Y.  Cheng  (1966). 

A  simplified  version  of  this  method  (known  as  the 
deghost  method)  which  neglects  noise  has  also  been  devel¬ 
opped  and  is  described  below.  Both  methods  have  been  com¬ 
pare  with  two  types  of  cross-correlation  analysis. 

Objective  5:  To  construct  a  two-dimensional  scale- 
model  of  the  earth's  crust,  to  examine  the  effects  of 
refraction  of  seismic  waves  along  the  Mohorovicic  dis¬ 
continuity  in  this  model  and  to  compare  them  with  results 
in  the  real  earth.  Technical  Report  Nr.  7  by  D,  E.  Siskind 
(1966)  describes  the  principal  experiment  in  this  series. 
Other  experiments  related  to  this  objective  are  described 
in  Technical  Report  Nr.  1  by  Y.  Nakamura  (1963)  and  in 
Scientific  Report  Nr.  1  (Nakamura,  1964),  Nr.  4  (Howell, 

1967)  and  Nr.  5  (Howell  and  Baybrook,  1967). 

Objective  6:  To  prepare  reports  covering  all  work 
done.  A  complete  list  of  all  reports  prepared,  including 
verbal  reports  before  societies,  published  reports  to  which 
this  contract  contributed  partial  support  and  publications 
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planned  which  will  be  completed  after  the  close  of  this 
contract  are  included  in  Table  I. 

■H*  Summary  of  Conclusions  arrived  at  in  Previously 
Issued  Reports. 

A.  Model  materials. 

In  the  early  experiments  on  this  project,  scale  models 
were  made  from  commercial  plastics.  These  materials  were 
found  to  absorb  elastic  waves  much  more  rapidly  than  metals. 
The  rates  of  absorption  were  measured  in  polystyrene  and  in 
Plexiglas.  Absorption  in  these  materials  was  found  to  be 

frequency  dependent.  The  absorption  coefficient  was  found 
to  be 

a  =  (0.054  log  F  -  0. 08o)/cm  (l) 

in  Plexiglas  and 

o  =  7-3  x  10  ^F^’^/cm  (2) 

in  polystyrene  for  30-180  kc,  where  F  is  frequency  and  a  is 
defined  by  the  equation  for  amplitude 

A  =  AQe-aA  (3) 

where  A  is  distance.  To  avoid  the  effects  of  absorption, 
only  metal  models  have  been  used  for  the  last  two  years. 
Absorption  is  much  lower  in  metals.  The  variation  with 

frequency  was  too  small  to  measure  in  the  bi-metallic  sheets 
used  in  this  research. 
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B.  Effect  on  the  frequency  spectrum  of  refraction 
along  boundaries. 

Three  types  of  boundary  were  investigated.  Nakamura 
(1963B)  studied  the  effect  of  a  transition  layer  in  which 
v/elocity  increased  smoothly  from  the  low-velocity  material 
to  the  high-velocity  material.  Seismic  pulses  refracted 
along  such  a  boundary  were  attenuated  about  12  db/octave 
above  a  frequency  given  by  1.7  times  the  velocity  gradient 
in  the  model.  Using  field  data,  Nakamura  and  Howell  (1964) 
and  Howell  (I966A)  found  a  similar  frequency  cut-off  effect 
corresponding  to  a  transition-layer  of  not  over  0.5  km 
thickness  under  the  Gulf  of  Maine  and  Lake  Superior. 

j.his  does  not  mean  that  such  a  transition  layer  exists 
in  these  places.  Siskind  (1966)  found  that  a  thin  inter¬ 
mediate-velocity  layer  could  have  an  effect  on  the  spectrum 
which  could  not  be  distinguished  from  a  Nakamura-tyoe 
transition  layer. 

Howell  and  Eaybrook  (1966)  investigated  an  irregular 
(saw-toothed)  boundary  (Fig.  2),  and  found  that  it  also  pro¬ 
duced  a  decrease  in  energy  at  high  frequencies  (Fig.  3).  It 
thus  appears  that  a  wide  range  of  structures  can  produce  the 
observed  effect,  which  must  consequently  be  expected  from 
many  if  not  nearly  all  real  crustal  structures.  Critical 
refraction  thru  a  layered  crust  or  along  the  Mohorovicic  dis¬ 
continuity  will,  in  this  fashion,  tend  to  limit  seismic 
signals  to  frequencies  below  6-10  cps.  The  peak  energy  in 
the  spectrum  tends  to  be  below  this.  This  filtering  effect 

/ 


will  have  a  tendency  to  obscure  any  differences  in  the 
spectrum  of  earthquakes  and  blasts,  even  if  such  differ¬ 
ences  existed  at  the  source. 

C.  Depth  of  focus  from  simple  spectral  analysis. 

Throughout  this  research  the  basic  principle  used  in 
an  attempt  to  measure  depth  of  focus  has  been  the  concept 
that  there  is  a  surface-reflected  pulse,  pP,  arriving  after 
an  unreflected  first  arrival,  P,  and  resembling  it  in  wave¬ 
form.  In  the  frequency  domain,  the  presence  of  these  two 
pulses  will  produce  maxima  and  minima,  the  positions  of  which 
are  determined  by  the  pP-P  delay  time.  Rothman  (1964)  sought 
to  measure  this  delay  by  examining  the  interference  patterns 
which  resulted  from  the  overlapping  of  these  two  pulses.  He 
found  maxima  and  minima  in  the  frequency  spectrum  at  the 
expected  places,  as  did  Siskind  (1966)  in  a  more  limited 
test.  Rothman  showed,  however,  that  even  in  the  case  of 
very  simple  layering,  there  were  many  superposed  interfer¬ 
ence  patterns  resulting  from  the  presence  not  only  of  P  and 
pP,  but  also  from  P  and  other  pulses.  Because  every 
passage  thru  an  Interface  potentially  produces  up  to  four 
pulses,  the  pattern  of  peaks  and  troughs  in  the  spectrum 
of  any  seismogram  from  a  real  structure  is  complicated. 
Rothman  showed  that  there  is  a  spectral  effect  which  marks 
the  depth  of  focus,  but  that  it  was  not  an  obvious  one 
which  could  be  measured  dependably  in  any  but  special, 
simple  cases.  Similar  results  were  obtained  by  autocorrela¬ 
tion  analysis. 
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Leblanc  (1966)  studied  the  complexities  of  the  spectrum 
further  ard  showed  that  even  from  one  earthquake  the 
spectrum  varied  greatly  from  station  to  station  (Figure  4). 
He  demonstrated  that  the  detailed  structure  of  the  spectrum 
can  be  fully  accounted  for  by  the  effects  of  reflection  and 
refraction  as  a  seismic  pulse  passes  thru  a  layered  crust. 

He  showed  that  this  effect  on  the  oulse,  called  the  transfer 
function  of  the  crust,  could  be  calculated  for  any  segment 
of  seismogram.  He  developped  a  method  of  choosing  between 
different  transfer'  functions  to  decide  which  of  two  or  more 
was  the  more  likely  cause  of  an  observed  spectrum  represent¬ 
ing  any  given  seismogram. 

If  the  structure  of  the  earth  ^s  accurately  known, 
Leblanc's  method  will  allow  the  spectrum  to  be  adjusted 
(thru  inverse  filtering)  until  only  the  effect  of  depth  of 
focus  remains.  However,  the  spectrum  is  so  sensitive  to 
small  changes  in  layer  thicknesses  and  physical  properties 
that  more  research  will  have  to  be  done  before  the  effect 
of  the  crust  at  any  station  can  be  removed  well  enough  to 
try  Rothman's  method  of  finding  focal  depth. 

D.  Watson's  focal-depth  method. 

Since  every  seismic  event  except  one  right  at  the 
surface  must  produce  both  a  P  and  a  pP  pulse  (except  for 
the  special  cases  of  zero  P  reflection  from  the  surface  and 
nulls  in  the  radiation  pattern  at  the  source)  the  seismogram 
normally  consists  of  several  trains  of  pulses:  l)  the  P 
pulse  travelling  from  the  source  to  the  recorder,  followed 
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by  all  the  additional  pulses  generated  by  reflection  and 
refraction  at  intermediate  boundaries;  2)  the  surface 
reflected  pP  pulse  and  its  train  of  additional  pulses; 

3)  the  surface  reflected  sP  pulse;  4)  the  S  pulse;  etc. 

If  only  that  portion  of  the  seismogram  before  the  arrival 
of  sP  is  considered,  this  portion  may  be  treated  as  the  s am 
of  rwo  trains  of  pulses:  P(t)  and  pP(t).  Where  the  upward- 
radiated  P  pulse  has  the  same  waveform  as  the  downward- 
radiated  P  pulse,  where  there  are  no  large  additional  pulses 
generated  by  Interfaces  between  the  focus  and  the  surface, 
and  where  the  angle  of  incidence  at  the  surface  is  small  so 
that  pP  and  P  travel  nearly  the  same  path  beneath  the  surface, 
then  the  seismogram  can  be  considered  to  consist  of  only  the 
sum  of  I (t),  pP(t)  and  any  noise  which  is  present.  Since 
pP(t)  is  postulated  to  be  of  the  same  waveform  as  P(t)  but 
delayed  behind  it  by  the  time.  At,  to  rise  tc  the  surface 
and  descend  again  to  the  focal  depth, 

pP(t)  =  RP(t-At)  (4) 

where  R  is  the  ratio  of  amplitudes  of  p?(t)  to  P(t).  In 
this  case,  the  seismogram  will  be 

S(t)  =  P(t)  +  RP(t-At)  +  N(t)  (5) 

where  N(t)  is  noise  plus  the  effects  of  imperfections  in  any 
of  the  approximations  assumed  above.  Knowing  S(t),  R  and 
At,  P(t)  can  be  approximated  for  any  noise  level. 

An  approximation  of  P(t)  can  be  found  for  any  S(t),  R 
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and  At.  If  R  and  At  are  known  to  lie  within  some  range  of 
values,,  these  can  be  sampled,  and  P(c)'s  found  for  all 
desired  values  of  R  and  At  for  any  given  S(t).  Merdler 
(  964)  experimented  with  different  criteria  for  picking  that 
P(t,R,At)  which  was  the  best  approximation  to  the  original 
pulse  radiated  from  the  focus  of  the  seismic  event.  He 
postulated,  on  empirical  evidence,  that,  for  the  correct 
values  of  At  and  R,  the  value  of  the  quantity  would  be 
larger  than  for  any  other  combination  of  At  and  R,  where 


2  (|p(t)|  -  U)2 

"pu _ 

"(2  t)-(ZP2(t)) 

T  ,  T 


(6) 


where  T  is  the  sequence  of  time  intervals  when  P(t)  exceeds 
P  ^ 

the  noise  level  U,  and  T  is  the  total  length  of  P(t)  tested. 

C_  is  called  concentration  ratio.  C„  was  considered  to  be  a 
n  n 

measure  of  the  relative  simplification  of  P(t)  as  compared 
to  S( t ) . 

Lin  ( 1966)  and  Cheng  (1966)  tested  the  Watson-Merdler 
method  on  fourteen  nuclear-test-blast  seismograms  and  on 
fourteen  earthquake  seismograms,  respectively.  They  found 
that  they  had  to  add  to  Merdler1 s  criterion  the  requirement 
that  the  numerator  in  equation  (4)  was  less  than  the  corre¬ 
sponding  summation  for  the  original  seismogram: 

2 
T 


(  |p(t)  I  -  U)2<  £  (  |S(t)  |  -  U): 

T 

pu 


su 


(7) 
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They  also  found  no  evidence  that  the  largest  C_  necessarily 
occurred  for  the  case  sought.  By  examining  all  P(t)'s  for 
combinations  of  R  u.nd  At  where  CR  was  greater  than  the  con¬ 
centration  ratio  for  the  original  seismogram,  S(t),  from 
which  P(t)  was  derived,  they  discovered  that  there  seemed 
to  be  more  than  one  P(t)  which  was  a  simplified  version  of 
S(t).  Each  was  suspected  of  being  a  case  where  the  inverse 
filtering  used  to  find  P(t)  in  the  Watson- Merdler  procedure 
had  removed  seme  real  pulse  or  pulses.  Sometimes  the  removed 
pulses  could  be  identified  as  pP,  sP  or  the  core  reflection, 
PcP.  At  other  times  they  could  not  be  identified. 

It  was  also  found  on  checking  that  the  effect  of  the 
layering  of  the  crust  was  to  produce  many  pulses  in  the 
seismogram.  It  was  concluded  that  the  most  likely  explana¬ 
tion  for  the  many  combinations  of  R  and  At  which  produced 
reasonably  simplified  p(t)*s  was  that  each  corresponded  to 
a  case  where  one  of  these  crustal  internal  reflections  was 
removed.  Which  P(t)  was  pP  was  nor  obvious. 

The  Watson-Meraler  method  succeeded  in  recognizing 
pulses,  but  it  could  not  identify  them  (Figs.  5  and  6). 

Even  in  the  case  of  a  shallow  blast,  there  were  usually 
several  P(t)’s  recognized.  Hence,  the  method,  in  its  ele¬ 
mentary  form,  failed  as  a  means  of  distinguishing  between 
blasts  and  seismograms.  Using  many  earthquake  seismograms, 
however,  there  is  a  good  possibility  that  pP  and  other 
pulses  can  be  recognized  by  their  tendency  to  occur  as 
solutions  regularly  for  a  large  fraction  of  the  records. 
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IV.  Comparison  of  the  Watson-Merdler  method  with 
alternative  approaches. 


The  Watson-Merdler  method  was  able  to  find  values  of 
time  delay.  At,  and  pulse-strength  ratio,  R,  for  which  the 
seismogram,  S(t),  seemed  to  consist  approximately  of  two 
pulses:  a  P(t),  which  was  simpler  than  the  original  seismo¬ 
gram  plus  an  image  of  itself,  RP(t-At).  Merdler's  measure 
of  simplicity,  CR,  was  largely  empirical.  In  the  course  of 
Lin's  and  Cheng's  studies,  it  became  increasingly  apparent 
that  the  use  of  CR  as  a  means  of  identifying  the  cases  where 
inverse  filtering  was  removing  a  real  pulse  was  the  weakest 
part  of  the  procedure.  C„  consists  of  three  parts: 

n 

Up  -.2  (|p(t)  I  -  u)2  (8) 

p  |p|>U  1 


where  U  is  the  noise  amplitude 


Lp  =  ZPc(t) 


T  =  .Z,  t 

p  |p|>u 


Using  this  notation,  equation  (4)  can  be  rewritten 


cR(p)  =  A- 

p  p 


(9) 

(10) 


(11) 


There  are  corresponding  quantities  for  S(t) 


U  =  2  (  |S(  t)  -  U): 

|s  |>u 


L  =  ZS  (t) 


(12) 


(13) 


-  12 


O' 

~s 


CR(S)  = 
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T  -L 


(14) 

(15) 


To  identify  cases  where  inverse  filtering  simplified  the 
seismogram,  Lin  (1966)  and  Cheng  (1966)  used  equation  7  and 
the  following  equation  in  place  of  the  condition  that  as 
given  by  equation  6  be  a  maximum 

CR(P)>  CR(S)  (16) 


The  second  test  was  necessary  because  in  certain  cases, 
inverse  filtering  produced  a  P(t)  whose  amplitude  was  much 
greater  than  that  of  S(t).  It  did  not  seem  reasonable  that 
P(t)  should  represent  more  energy  than  the  S(t)  of  which  it 
was  presumed  to  be  a  part. 

Peaks  in  the  Watson-Merdler  CR  occur  largely  when  Tp  is 

significantly  less  than  T  .  The  two  criterion  T  <T  and 

s  p  s 

can  be  combined  to  give  a  test  criterion 

p  s 


>1 


(17) 


This  test  was  proposed  late  in  the  research  as  a  sub¬ 
stitute  for  the  Watson-Merdler  CR  test.  It  has  not  been 
properly  tested  to  see  if  it  is  a  better  test.  It  was  used 
in  place  of  Cr  in  testing  the  simplified  ^rsion  of  the 
Watson-Merdler  procedure  described  below. 

The  Watson-Merdler  method  finds  p(t)'s  which  meet 
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certain  criteria  within  a  designated  noise  level,  U.  If  it 
is  assumed  that  no  noise  is  present;,  ohen  a  P(t)  can  still 
always  be  found  for  any  S(t)  for  any  combination  of  At  and 
R  for  a  finite  time  after  the  start  of  S(t)  provided  S(t) 
is  zero  before  time  t=0.  In  this  case  equation  (3)  becomes 

S(t)  =  P(t)  +  RP(t-At)  (18) 

This  can  be  rearranged,  giving 

P(t)  =  S(t)  -  RP(t-At)  (19) 


It  is  simple  to  write  a  computer  program  which  will  find 
P(t)  from  S(t),  R  and  At.  This  was  done.  This  program  is 
called  the  "deghost"  program.  This  is  possible  because 
P(t)  =  S(t)  up  to  t=A t .  For  every  t  between  At  and  2At, 

P(t)  can,  therefore,  be  calculated  from  already  known  values 
of  S(t).  After  calculating  P(t)  for  the  range  Atto2Ao-,  the 
calculation  can  be  extended  to  3At,  and  so  forth. 

Such  a  program  has  one  major  weakness:  the  errors  due 
to  noise  can  be  cumulative,  and  the  pulse  can  increase  in 
average  amplitude  with  time.  An  S(t)  which  gives  an  unstable 
P(t)  is  the  continuous  signal 


=  0 

S( t)  =  ln( 1  +  at ) 


O  ^  t 

0<t£At 


=  ln( 


1  +  a  t 


1+  at  -a At 


■) 


At£t 


(20) 


Solving  ( 18)  for  P(  t)  when  R=-l  and  At£t<2At 


-  14  - 


P(t)  =  S(t)  +  P(t-At) 

“  ‘"(kUti  +  ln(l-wt-aAt) 

=  ln(l-fat)  (21) 

so  that 

P(t)  =  0  t<0 

=  ln(l+at)  t<0  (22) 

which  rises  continuously  toward  infinity  as  At  increases. 
Provided  aAt,  is  less  than  unity,  S(t)  increases  from  zero 
at  t*0  to  ln(l+at)  at  t=At,  and  thereafter  decreases  to  zero 
as  t  approaches  infinity.  Any  noise  pulse  which  resembles 
this  will  produce  a  component  of  P(t)  which  steadily 
increases  in  size. 

The  deghost  program  was  tried  on  the  same  seismograms 
used  to  test  the  Watson-Merdler  program.  A  similar  procedure 
has  been  used  by  Sarmah  and  Berg  (1966),  who  credit  it  to 
Carpenter  (1964).  Computations  by  the  deghost  method  are 
accomplished  much  faster  than  by  the  Watson-Merdler  method, 
using  less  than  one  quarter  of  the  computer  time. 

Tables  II  and  III  compare  the  cases  of  greatest  simpli¬ 
fication  found  by  the  Watson-Merdler  method  with  those 
found  by  the  deghost  method.  CK  was  used  in  place  of  CR  as 
a  criterion  for  simplification  in  the  deghost  program.  For 
the  fourteen  earthquake  seismograms,  the  deghost  method 
found  that  a  simplification,  as  indicated  by  a  peak  in  CR, 
had  occurred  for  At's  within  0.15  seconds  of  the  greatest 
simplification  found  by  the  Watson-Merdler  method  in  only 
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seven  nut  of  the  17  cases  of  strongest  Watson-Merdler  simpli¬ 
fication,  testing  over  the  same  range  in  both  cases.  (The 
deghost  program  was  run  over  a  wider  range  of  time  delays 
than  the  Watson-Merdler  method.  The  largest  deghost  CR's 
occurred  in  some  cases  for  delays  shorter  than  those  tested 
by  the  Watson-Merdler  method.)  Since  the  P(t)'s  found  by 
the  deghost  method  have  not  been  systematically  studied,  it 
is  unknown  which  of  the  two  methods  is  the  better  at  picking 
the  most  obvious  cases  of  simplification.  Nor  is  it  known 
to  what  degree  the  differences  in  the  results  are  due  to 
differences  between  the  CR  and  the  CR  test  criteria  and  to 
what  degree  they  reflect  the  failure  of  the  deghost  program 

to  allow  for  the  effects  of  noise. 

In  the  case  of  the  fourteen  nuclear  blast  seismograms,, 
the  deghost  method  indicated  simplification  in  ten  out  of 
twelve  of  the  cases  of  greatest  CR  picked  by  the  Watson- 
Merdler  method;  and  the  Watson-Merdler  method  found  simpli¬ 
fications  in  nine  out  of  fourteen  of  the  cases  of  greatest 
CTJ  found  by  the  deghost  method  (Table  III). 

The  first  parts  of  the  fourteen  blast  and  fourteen 
earthquake  seismograms  were  also  cross -correlated  with  the 
whole  seismograms,  and  the  cases  of  strongest  correlation 
are  reported  in  Tables  II  and  III.  Strong  correlations  with 
the  greatest  simplifications  picked  by  the  Watson-Merdler 
method  occurred  in  10  out  of  28  cases,  and  with  the  deghost 
largest  CH's  in  15  out  of  23  cases. 

One  other  test  was  made.  Define  a  segment  of  a  digitized 

seismogram  from  the  A’th  to  the  B'th  point  as  S^_g(-)* 


. .  'Il’l' .  . . .  . ..,.UMt.WMl.,»n)|„HpulmHUtMil1^1 
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Define  another  segment  of  equal  length  starting  with  the 
C'th  point  as  SQ_D(t),  where  D-C=B-A.  Now  define  Qc(t)  as 

V*)  -  Sc-D(t)  -  “a-b^)  (23) 

where  R  is  a  multiple  of  0.1.  The  "search"  program  com¬ 
pares  Sc_D(t)  with  SA_B(t)  as  a  function  of  the  variable,  C, 
by  finding  the  largest  value  of  +R  or  -R  where  Q  (t)  is 
simpler  than  SA_R(t)  using  condition  17.  This  method  expeot- 
ably  gives  very  nearly  the  same  results  as  cross-correlation. 

These  tests  are  not  conclusive.  Examination  of  the 
results  has,  however,  led  to  the  suspicion  that  the  basic 
trouble  may  not  lie  in  the  Watson-Merdler  or  deghost  methods 
themselves,  but  in  the  nature  of  the  seismograms.  The  seis¬ 
mograms  used  were  all  recorded  with  short-period  Benioff  or 
similar  high-natural -frequency  seismometers.  The  displace¬ 
ment  sensitivity  of  the  Benioff  instrument  falls  off  rapidly 
below  its  peak  sensitivity  at  1.67  cycles  per  second,  and  is 
about  24  db  down  by  0.4  cps .  The  energy  in  the  earthquake 
studied  by  Leblanc  (Technical  Report  Nr.  4,  1966)  usually 
had  its  peak  at  or  below  0.4  cps.  The  effect  of  this  is  to 
make  the  energy  in  the  seismograms  sharply  peaked  in  fre¬ 
quency  near  1  cps.  This  is  born  out  by  a  tendency  for  peaks 
in  the  cross-correlation  function  to  occur  at  regular  one- 
second  intervals.  This  is  true  also  of  the  deghost  solutions. 
This  was  not  so  apparent  in  the  case  ox’  the  Watson-Merdler 
method,  but  whether  this  is  because  of  a  greater  sensitivity 
of  that  method,  or  because  the  CR  test  was  poorer  than  the 
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CH  test  for  picking  simplifications,  or  for  some  other 
reaso  .  is  unknown.  One  thing  is  apparent,  however:  the 
seismic  data  analyzed  should  include  lower  frequencies 
than  the  predominant  range  in  the  seismograms  tested  here. 

Another  reason  why  the  method  may  have  failed  to  pick 
pP  consistently  in  seismograms  is  that  the  pulse  shape  of 
pP  is  not  always  the  same  as  that  of  p.  To  test  this, 
several  earthquakes  of  presumably  known  focal  depth  have 
been  examined.  Figure  7  is  the  first  part  of  the  State 
College,  Pa.,  seismogram  of  the  11  January  1965  Alaskan  earth¬ 
quake,  reported  by  the  Coast  and  Geodetic  Survey  to  be  at  59 
km  focal  depth.  This  focal  depth  corresponds  to  a  pp-p  delay 
of  about  13.5  sec.  There  is  no  clear  pulse  at  this  time  in 
the  seismogram,  but  there  is  a  strong  pulse  beginning  at  7.6 
sec  (corresponding  to  31  km  focal  depth).  The  deghost  pro¬ 
gram  finds  this  correlation  clearly.  It  also  finds  a  weak 
correlation  at  9-15  sec  (38  km).  It  finds  no  simplification 

for  R=-0. 5  for  delays  between  9-15  and  25  sec,  the  largest 
delay  tested. 

Figure  8  compares  the  spectra  of  the  Initial  P  and  the 
pulse  at  7.6  sec.  Although  these  are  similar,  they  are  far 
from  identical.  Leblanc's  study  (Technical  Report  Nr.  4) 
of  tne  crustal  transfer  function  shows  that  the  crust  in 
the  epicentral  region  is  capable  of  introducing  enough 
changes  into  the  pulse  shape  of  pP  to  cause  it  to  differ 
substantially  from  P  in  the  frequency  band  used  here. 

There  is  also  the  possibility  that  the  ghost  intensity 

% 
fa  - 
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coefficient,  R,  Is  sometimes  positive  for*  pP  for  earthquakes. 
The  first  motion  from  earthquakes  can  be  a  compression  In 
some  azimuths,  a  dilatation  in  others.  Neither  Lin  nor 
Cheng  tested  systematically  for  positive  values  of  R. 

Tests  using  the  deghost  program  found  approximately  as  many 
cases  of  CH>  1  for  positive  R  as  for  negative  R  in  a  few 
trial  runs.  Large  positive  cross-correlations  also  occurred 
about  as  often  as  negative  ones. 

V*  Conclusions  and  Recommendations 

This  research  has  demonstrated  that  the  Watsen-Merdler 
and  the  deghost  methods  can  find  pulses  in  a  seismogram,  but 
that  in  the  frequency  range  recorded  well  by  short-period 
seismometers,  more  pulses  are  found  than  can  be  easily  inter¬ 
preted.  The  research  suggests  but  has  not  proven  that  pulse 
interference  due  to  the  crustal  transfer  function  is  a  major 
disturbing  factor  and  may  be  the  principal  cause  of  this 
difficulty.  As  a  result  of  this,  these  methods  have  failed 

to  provide  an  easily  used  means  of  distinguishing  earthquakes 
from  blasts. 

One  of  the  major  needs  in  seismology  is  a  means  of 
removing  the  effects  of  the  crustal  transfer  function  from 
seismograms.  It  appears  likely  that  it  is  the  complexities 
introduced  by  crustal  layering  that  have  prevented  these 
methods  from  separating  blasts  from  earthquakes  and  have 
Introduced  uncertainties  into  the  depth-of-focus  measure¬ 
ments.  Leblanc  (Technical  Report  Nr.  4,  1966)  has  suggested 
an  approach  by  which  this  problem  may  be  solved.  When  means 
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of  removing  the  crustal  effect  from  the  seismogram  at  a 
recording  station  have  been  developped,  this  whole  depth- 
of-focus  experiment  should  be  repeated. 

This  research  has  shown  also  that  seismic  energy  falls 
off  rapidly  above  0.4  cps.  Any  further  work  done  along  the 
lines  tried  here  should  use  data  covering  a  broader  range 
of  frequencies  at  the  low-frequency  end. 

The  methods  developped  here  do  show  promise  of  provid¬ 
ing  a  means  of  selecting  the  pP-P  time  interval  more  accurately 
than  by  simple  visual  inspection  of  the  seismogram.  This 
possibility  should  be  studied  in  greater  detail  on  longer- 
period  seismograms.  Long  period  seismograms  are  advantageous 
also  because  pulse  interference  is  less  prominent  than,  in 
short-period  seismograms. 

The  Watson-Merdler  and  deghost  methods  should  be  com¬ 
pared  more  thoroughly  to  find  if  the  complicated  Watson- 
Merdler  method  is  significantly  more  sensitive  to  ghosts  in 
a  seismogram  than  is  the  deghost  method,  and  to  find  out  if 
either  is  better  at  picking  ghosts  than  cross-correlation  or 
autocorrelation. 
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„°  parts  of  the  scale-model  labor 
Upper  view:  a  model  in  place  ready 
pulses.  Lower  view:  apparatus  for 
earthquake  seismograms. 


Fig.  2.  Model  used  to  investigate  irregular  Mohorovicic 
boundary. 
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FREQUENCY  (KC) 

Fig.  3.  Spectra  of  pulses  critically  refracted  at  two 
distances  along  irregular  boundary  compared  to 
spectrum  of  same  pulse  before  refraction  (at 
14  cm) . 


5  10  15  20  a«c 

Pig.  5-  Original  3(t)  (solid  line)  and  deconvolved 
P(t)  (dashed  line)  earthquake  seismograms 
whi<_':  show  simplification  of  P(t)  as  compared 
to  S^O*  Ecuador  Earthquake  of  10  May  1963. 


Fig.  6.  Original  S(t)  (solid  line)  and  deconvolved 
P(t)  (dashed  line)  blast  seismograms  which 
show  pulse  removal  by  inverse  filtering. 
Salmon  blast  at  Mould  bay,  Canada. 
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